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ABSTRACT 

We report the detection of extended emission around the anomalous X-ray pulsar AXP 1E1547.0- 
5408 using archival data of the Chandra X-ray satellite. The extended emission consists of an inner 
part, with an extent of 45" and an outer part with an outer radius of 2.9', which coincides with a 
supernova remnant shell previously detected in the radio. We argue that the extended emission in 
the inner part is the result of a pulsar wind nebula, which would be the first detected pulsar wind 
nebula around a magnetar candidate. Its ratio of X-ray luminosity versus pulsar spin-down power is 
comparable to that of other young pulsar wind nebulae, but its X-ray spectrum is steeper than most 
pulsar wind nebulae. We discuss the importance of this source in the context of magnetar evolution. 
Subject headings: starsmeutron - starsrmagnetic fields - ISM:supernova remnants - 
(stars:)pulsars:individual (AXP 1E1547.0-5408) 



1. INTRODUCTION 

Among the various types of neutron stars, soft gamma- 
ray repeaters (SGRs) and anomalous X-ray pulsars 
(AXP) stand out by their extreme bursting behavior, 
high X-ray luminosity, and long rotation periods of 2- 
12 s (see Woods & Thompson 2004; Mereghetti 2008, 
for reviews). They are now widely recognized as magne- 
tars, neutron stars with extreme surface magnetic fields 
of 10^'' — 10^^ G, well above the quantum critical mag- 
netic field of Bqed = 4.4 x 10^^ G. The X-ray luminosity 
and bursts of magnetars are thought to be powered by 
the decaying magnetic field. Despite their extreme en- 
ergy output, the particle acceleration properties of mag- 
netars are not well known. 

Here we report the detection of a pulsar wind neb- 
ula (PWN) around AXP IE 1547.0-5408. We show that 
the nebula's brightness is comparable to that of normal 
pulsars with the same rotational energy loss rate. This 
suggests that the nebula is powered by rotational energy 
loss, rather than by the magnetar's magnetic field. 

1.1. IE I547.O-54O8 

IE 1547.0-5408 (Lamb & Markert 1981) is an X-ray 
source that was only recently recognized to belong to 
the class of magnetars (Gelfand & Gaensler 2007), a 
name we adopt here to refer to both SGRs and AXPs. 
With its rotation period of P = 2.1 s (Camilo et al. 
2007), it is the most rapidly spinning magnetar known. 
Its period derivative is P = 2.3 x lO"^'^, which com- 
bined with the period, implies a surface magnetic field 

of Bd^pole = 3.2 X 10^^ ^/pP = 2.2 X 10^^ G, and a ro- 
tational energy loss rate of -E = 1.0 x 10'^^ ergs~^, the 
highest among AXPs and SGRs. IE 1547.0-5408 is also 
a transient radio source (Camilo et al. 2007), a charac- 
teristic shared with only one other magnetar, AXP XTE 
J1810-197 (Camilo et al. 2006). Very recently IE 1547.0- 
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5408 went through a period of major bursting activity 
(Krimm et al. 2008b, a; Gronwall et al. 2008), confirming 
its status as a magnetar. 

IE 1547.0-5408 is located inside a supernova remnant 
(SNR) shell, G327.24-0.13, detected with the Molongo 
radio observatory (Gelfand & Gaensler 2007). Although 
the spin-down rates of most AXPs/SGRs imply that they 
are young sources (typically < 10^ yr), only a few of them 
are associated with SNRs (Gaensler 2004). These SNRs 
may hold important clues about the supernovae that cre- 
ated the magnetars. In particular, the energy of the rem- 
nants indicates that the supernovae creating magnetars 
are not more powerful than ordinary supernovae (Vink 
& Kuiper 2006). This poses a challenge to the scenario 
in which magnetars are born with extremely high initial 
spin-periods (Duncan & Thompson 1992), in which case 
10^^ erg rotational energy is likely to be transfered 
to the supernova ejecta. Nevertheless, this scenario can- 
not be completely excluded, because the distances used 
by Vink & Kuiper (2006) have been disputed (Durant & 
van Kerkwijk 2006). 

A major obstacle in determining the initial spin- 
periods of AXPs/SGRs is that they spin-down rapidly 
and often irregularly (e.g. Gavriil et al. 2004; Espos- 
ito et al. 2009). It is, therefore, difhcult to determine 
whether they were born with periods close to their cur- 
rent spin-periods, or with periods comparable to, or even 
shorter than ordinary pulsars. Provided that the nebula 
around IE 1547.0-5408 is indeed powered by rotational 
energy loss, its properties may help to clarify its rota- 
tional history, as PWNe contain most of the initial ro- 
tational energy lost over the lifetime of the pulsar (e.g. 
Pacini & Salvati 1973). 

In the next section we report on the details of the 
detection of extended emission around IE 1547.0-5408, 
both based on imaging (§ 2.2), and on spectral analysis 
(§ 2.3). In § 2.4 we arguethat the extended source is most 
likely a PWN, and not a dust scattering halo. In the final 
section we discuss the implications of the presence of a 
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Fig. 1. — Left: A Chandra ACIS-I image of the diffuse emission surrounding IE 1547.0-5408 in the energy range 1.5-7 keV. Several point 
sources, including IE 1547.0-5408, were removed, and the image was then smoothed with a kernel of cr = 11.8" (24 pixels). Right: The 
radial brightness profile in the 1.5-7 keV band (black squares) and radio (dotted line, based on the SUMMS survey (Bock et al. 1999)). We 
have normalized the radio profile such that the central enhancement coincides with extended emission in X-rays. The X-ray background 
intensity is indicated by the dashed line, and the spectral extraction regions for the PWN, shell and background (bgd) by two-sided arrows. 
With cyan triangles we also indicate the profile in the 3-7 keV range. 



PWN around IE 1547.0-5408. 

2. DATA ANALYSIS AND RESULTS 

2.1. Data selection 

In order to find a potential X-ray counterpart to the 
radio shell of G327.24-0.13, we analyzed archival Chan- 
dra and XMM-Newton data. During this search we found 
a faint, but distinct, extended X-ray source, centered on 
the magnetar. Not all archival data were suitable for our 
analysis, since several Chandra observations were taken 
in continuous clocking mode, which does not allow for 
imaging. This leaves only the 9.6 ks long observation 
on July 1, 2006, with the ACIS-I detector, as a suitable 
observation for extended source analysis. Several XMM- 
Newton observations are present in the archive, but the 
point spread function (PSF) of XMM-Newton is much 
broader ('^ 6" FWHM, with extended broad wings) than 
that of Chandra (- 0.5" FWHM). Therefore, the Chan- 
dra data, although poor in statistics, provide a more con- 
vincing case for the presence of an extended source. How- 
ever, we note that the XMM-Newton data are consistent 
with the results presented here. 

2.2. Imaging 

We performed the analysis of the Chandra data using 
the standard software (CIAO v4.1) to make exposure cor- 
rected images, and extract spectra of the magnetar, the 
extended emission, and the region of the radio shell that 
was found by Gelfand & Gaensler (2007). The extended 
source is weak, and can only be seen after significantly 
smoothing the images. However, the advantage of the 
small PSF of Chandra is that point sources can be ef- 
fectively removed before smoothing. Apart from X-ray 
emission within 4" of the position of IE 1547.0-5408 also 
other faint point sources in the field were removed, two 
of which were located in the region of the radio shell^. 



1 Their J2000 positions are a = Ib^^bV 
and a = 15''50™34.8^, 5 = -54°19'25.9" 



'07.7^5 = -54°19'26.0", 



The resulting 1.5-7 keV image is shown in Fig. 1 (left). 
We chose for this energy range, because below 1.5 keV 
there is hardly any emission to be expected due to to 
de large absorption column (c.f. Fig. 2 in Gelfand & 
Gaensler 2007). The image is noisy, but extended emis- 
sion centered on IE 1547.0-5408 can be seen. This source 
is surrounded by a faint shell. 

The radial X-ray brightness profile makes the struc- 
ture even more clear (Fig. 1 right). We include in this 
figure the radio intensity profile at 843 MHz as obtained 
from the Sydney University Molonglo Sky Survey (Bock 
et al. 1999; Gelfand & Gaensler 2007). We also show 
the profile using the hardest photons only, 3-7 keV, as 
this lowers the potential contribution of dust scattered 
X-ray emission from IE 1547.0-5408 (see § 2.4). The pro- 
file clearly shows the presence of an extended emission 
component that falls off with radius and disappears at 
r w 45". Further outward there is a small enhancement 
out to r « 150", which coincides with the radio shell. 
Also the radio map gives an hint of centrally enhanced 
emission with an approximate flux density of ^0.008 Jy 
at 843 MHz and a radial profile similar to that in X- 
rays, but with a different brightness ratio between shell 
and central emission. 

The radial brightness profiles allow us to make an es- 
timate of the statistical detection significances. Fitting 
the first three bins, labeled PWN, with a flat background 
model results in x^/d.o.f = 55.0/3, corresponding to a 
6.9(7 detection. In the 3-7 keV band this reduces to 
X^/d.o.f = 15.3/3, still a 3.2cr detection. For the ra- 
dio shell the numbers are x^/d.o.f = 58.6/7 (6(t) and 
xVd.o.f = 18.1/7 (2.5cr), respectively. 

2.3. Spectral analysis 

We extracted X-ray spectra of three regions, using a 
radius of r = 0.74" to isolate the spectrum of IE 1547.0- 
5408 (corresponding to 80% encircled energy); an annu- 
lus with inner radius of 4" and an outer radius of 45" for 
the extended source; and an annulus between 45" and 
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Fig. 2. — The Chandra X-ray spectra of the pulsar wind nebula, 
the supernova remnant shell and IE 1547.0-5408. The countrates 
of IE 1547.0-5408 and the pulsar wind nebula have been multiplied 
by the factors indicated for reasons of clarity. The rebinning was 
done for display purposes. No rebinning was necessary for the 
fitting itself as the statistic described in Cash (1979) was used. 



TABLE 1 

Spectral fit parameters^ 



Component 


Flux (2-10 keV) 


Spectral index 




(10~^^ erg s~^cm~-^) 






absorbed'' unabsorbed 




AXP 


2.4 ±0.3 3.8 ±0.4 


-3.8 ±0.1 


Extended source 


1.0 ±0.3 1.5 ±0.3 


-3.4 ±0.4 


Shell 


2.3 ±0.7 3.5 ±0.7 


-3.3 ±0.8 



''Based on joint fits of the spectra in the 1.5-7 keV range using 
the xpsec vll package, and employing the C-statistic, (Cash 1979). 
The overall statistic for the joint fit was C = 74.5 for 72 bins. 
Excluding the PWN model fit resulted in AC = 125 (Her), and 
excluding the Shell model resulted in AC = 110 (6a-). 

''The best fit absorption column was Nh = (2.75 ± 0.25) X 
10-^^ cm"-^, using the tbabs absorption model (Wilms et al. 2000). 

174" for the X-ray counterpart of the radio shell. The 
inner radius for the extended source guarantees that only 
~ 5% of the point source contributes to the emission of 
the extended source, due to the scattering wings of the 
Chandra PSF. ^ Given that the flux of IE 1547.0-5408 ap- 
pears to be about twice as large as of the extended source 
(Table 1), we expect a ~ 10% contamination of the ex- 
tended source from the AXP, due to the wings of the 
PSF. The background spectrum of the extended source 
and the shell was obtained from an annulus centered on 
the AXP with 180" <r < 272". The background spec- 
trum for the AXP was obtained from an annulus with 
4" < r < 20", i.e. we treated the extended region as a 
source of background for the point source. 

We fitted these data with power law spectra with a 
Galactic absorption component. For our analysis we 
also investigated an additional black-body component for 
the point source, and a non-equilibrium ionization thin 
plasma model for the radio shell. However, this did not 
improve the fits. Note that for the XMM-Newton data 

The Chandra Proposers' Observatory Guide Vll.O Fig. 4.20 
shows that 95% of the photons with energies of 1.5 keV fall within 
2" of the source position. 
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"Prediction based on Eq. 19 of Draine (2003). 
/3 = corresponds to uniformly distributed gas, /3 = 1 has 3/4th 
of the dust at a distance greater than 50% of the total distance 
(Draine 2003). 

including a black-body for the point source did improve 
the fit (Gelfand & Gaensler 2007), but the XMM-Newton 
data has better statistics concerning the spectrum of the 
point source. According to Gelfand & Gaensler (2007) 
the best spectrum of the AXP is a combination of a soft 
power law with index -3.7 plus a black-body component 
with kT = 0.43 keV. The Chandra data agrees with this 
in the sense that we also find a steep power law. The 
resulting spectra are shown in Fig. 2 and the best fit 
parameters are listed in the Table 1. 

2.4. A pulsar wind nebula or a dust scattering halo? 

The spatial morphology and the spectral characteris- 
tics of the extended source have all the characteristics of 
a PWN created by IE 1547.0-5408, but we need to ex- 
clude the possibility that the extended emission could be 
caused by a dust scattering halo (e.g. Predehl & Schmitt 
1995). Indeed, dust scattering is not absent as indi- 
cated by the recent detections of evolving dust scatter- 
ing haloes after some major outbursts of IE 1547.0-5408 
(Tiengo et al. 2009). However, in the situation investi- 
gated here a dust scattering halo is unlikely to dominate 
the emission above 3 keV, with a most likely fiux contri- 
bution to the emission within 45"of 6-11%. 

Our estimate is based on the relation between the dust 
scattering optical depth Tsca at 1 keV and the absorp- 
tion column as found by Predehl & Schmitt (1995). For 
iVn = 3 X 10^2 cm-2 (Table 1 and Gelfand & Gaensler 
2007) this gives Tgca — 1-5 at 1 keV, corresponding to 
Tsca = 0.17 at 3 keV, because the dust scattering cross 
section is proportional to E~^. The values of Tgca can 
be related to the fractional halo intensity, Ifmc (Predehl 
& Schmitt 1995). At 3 keV this gives Ifrac = 15%, but 
most of the halo fiux will be outside the spectral extrac- 
tion radius of 45". To estimate the fraction within 45", 
we use the semi-analytic approximation of Draine (2003), 
which employs a parameter (3 to account for gradients 
in the dust distribution; with f3 = ^l(l) correspond- 
ing to a dust distribution skewed in distance toward the 
observer (source), = to a uniform distribution. Ta- 
ble 2 lists the potential dust halo contributions to the 
extended source, which indicates a modest contribution 
of 11% even for /3 = 1. 

Draine (2003) discussed the possibility that the rela- 
tion iVn versus Tsca relation of Predehl & Schmitt (1995) 
underestimates Tsca, due to the neglect of small an- 
gle scattering. For that reason we also provide in Ta- 
ble 2 values for the expected dust halo contributions for 
Tsca = 3.0. In that case the dust halo becomes more sub- 
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Fig. 3. — The X-ray luminosity, Lx{2 — lOfcey), and the ro- 
tational energy loss, Erot of the magnetar IE 1547.0-5408 (filled 
triangle) compared to the unpulsed X-ray emission of other X-ray 
pulsars (Cheng et al. 2004). A distance of 9 kpc for IE 1547.0-5408 
is assumed (Camilo et al. 2007). The dashed lines indicate the 
standard deviation around the linear regression relation between 
Erot and Lx- Also indicated are the high magnetic field pulsars 
PSR J1846-0258 (Kargaltsev & Pavlov 2008; Gavriil et al. 2008) 
and PSR J1119-6127 (Safi-Harb & Kumar 2008) with filled circles, 
and RRAT J1819-1458 (Rea et al. 2009) with a star-like symbol. 

stantial, but at 3 keV, with 25% for P — 1, it still not 
dominant. Note that the model by Draine (2003) seems 
to overestimate dust halo scattering: Apart from its dis- 
crepancy with the results of Predehl & Schmitt (1995), 
also a recent Chandra observation of GX13-I-1 shows that 
this model gives a poor fit to the data and underestimates 
the absorption column (Smith 2008), and overestimates 
the fraction of small angle scattering. This suggest that 
the 25% contribution at 3 keV is too conservative, and an 
< 11% seems likely. Moreover, adopting Tgca = 3.0 pre- 
dicts a steepening of the halo spectrum within 45" with 
respect to the source spectrum of AF = —0.9, whereas 
the actual spectrum of the extended source is flatter than 
that of the AXP by AT = 0.4 ± 0.5 (a 2.7cr difference). 
Note that the expected steepening, is determined by the 
interplay between the energy dependence of the cross sec- 
tion, a oc and the energy dependence of the scat- 
tering angles. Still higher values of r^ca result in even 
stronger spectral steepening. 

3. DISCUSSION AND CONCLUSIONS 

We report here the presence of an extended source 
around IE 1547.0-5408, which is most likely a PWN. Not 
counting the spurious (Kaplan et al. 2002) association of 
a radio nebula with SGR 1806-20 (Kulkarni et al. 1994), 
this is the first evidence for a PWN around a magnetar. 
This is of particular interest, as young PWNe are formed 
during the whole lifetime of the pulsar, and thus hold a 
record of the energy released over the lifetime of the pul- 
sars (Pacini & Salvati 1973; van der Swaluw & Wu 2001; 
Chevalier 2004). This is important for magnetars, be- 
cause according to one theory magnetars are born with 
rapid initial spin periods (Pi < 3 ms) (Duncan & Thomp- 
son 1992, 1996), a theory that was recently challenged on 
observational grounds (Ferrario & Wickramasinghe 2006; 
Vink & Kuiper 2006) . But even if initial periods shorter 
than 3 ms seem unlikely, little is known about the initial 



spin-periods of magnetars. The presence of a PWN may 
therefore cast new light on this issue, as it may provide 
evidence that the initial spin period was much shorter 
than the present spin-period. 

In addition, the PWN around IE 1547.0-5408 may pro- 
vide new insight into the connection between AXP/SGRs 
and other young pulsars that have slightly lower mag- 
netic fields (McLaughlin et al. 2003; Kaspi & McLaughlin 
2005), and that do not have all the magnetar character- 
istics. A recent AXP-like burst from the rotation pow- 
ered pulsar PSR J1846-0258 already suggests that there 
is a continuum in properties from pulsars with average 
magnetic fields of 10^^ G to magnetars (Gavriil et al. 
2008). PSR J1846-0258 has an inferred surface magnetic 
field of 4.9 x 10^^ G and is surrounded by a PWN and 
a SNR sheU, Kes 75 (Kargaltsev & Pavlov 2008). An- 
other pulsar with a relatively high surface magnetic field 
{B = Ax 10" G), PSR J1119-6127, is also surrounded by 
a PWN, which appears jet-like in morphology (Safi-Harb 
& Kumar 2008). A third case of a high magnetic field 
pulsar with a PWN, RRAT J1819-1448, was recently re- 
ported by Rea et al. (2009). 

One problem in interpreting the presence of a PWN 
around a magnetar is that we do not a priori know 
whether the PWN is powered by rotational energy loss, 
as is the case for normal pulsars, or whether the neb- 
ula is ultimately powered by the magnetic field of the 
magnetar. For example it has been suggested that mag- 
netar magnetospheres are filled with electrons/positrons 
and/or ions (Thompson et al. 2002), a fraction of which 
may escape in the form of a relativistic particle wind 
(Harding 1996; Harding et al. 1999). 

A hint that the PWN around IE 1547.0-5408 may be 
powered by the pulsar spin-down is provived by com- 
paring the ratio of the X-ray luminosity over the spin- 
down power of the PWN. There is a correlation between 
the X-ray luminosity and spin-down power of normal 
young pulsar (Seward & Wang 1988; Verbunt et al. 1996; 

Cheng et al. 2004), with rjx — Lx/Erot in the range of 
lO"** — 10~^. Fig. 3 shows the correlation as taken from 
the sample of Cheng et al. (2004). It shows that the 
PWN around IE 1547.0-5408 is consistent with the gen- 
eral trend, with rj « 0.01, although perhaps with a some- 
what enhanced X-ray luminosity compared to pulsars 
with a similar rotational energy loss (at the 1.5(7 level). 
This is unlike the case of RRAT J1819-1448, which has 
rj = 0.2 (Rea et al. 2009). The value of 0.01 for 
IE 1547.0-5408 is consistent with the idea that its PWN 
is powered by the pulsar's spin-down. Interestingly, this 
can be taken as an additional argument against the al- 
ternative theory for the behavior of AXP/SGRs namely 
that they are not magnetars, but neutron stars that have 
a high spin-down rate due the propeller mechanism (e.g. 
Marsden et al. 2001). 

The idea that the PWN is powered by the pulsar spin- 
down is also credible given the fact that IE 1547.0-5408 
has the largest spin-down power among AXPs/SGRs. If 
the PWNe around magnetars are powered by magnetic 
activity, it is not clear why other AXPs/SGRs, some 
even more active than IE 1547.0-5408, do not show ev- 
idence for extended emission. On the other hand, it is 
not quite clear why the X-ray spectrum of the PWN is 
rather steep, F « —3.5, compared to PWNe around ordi- 
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nary pulsars, which have T w —2 (Kargaltscv & Pavlov 
2008). This may be taken as an hint for another origin, 
but it may also be explained by the much faster spin 
evolution of a magnetar compared to other young pul- 
sars. Using the terminology of Pacini & Salvati (1973), 
most young PWN are in phase 2 of their evolution, i.e. 
the pulsar has been a more or less constant source of 
relativistiv particles. In that case the spectral index be- 
yond the synchrotron break is steeper by only AP = 0.5 
compared to the radio spectral index. The PWN around 
IE 1547.0-5408 is most likely in phase 3: its spin-down 
luminosity has been rapidly decreasing during its short 
life {Tchar ^ 1400 yr). This means that the average age of 
the relativistic electron populations is skewed toward the 
total age of the PWN. However, this does not necessarily 
lead to a steeper spectrum (e.g. Pacini & Salvati 1973). 
This is an aspect of the nebula around IE 1547.0-5408 
that needs further clarification. 
Finally, we emphasize that our discovery is based on a 



Chandra observation of only 9.6 ks. A much deeper ob- 
servation is likely to reveal further details. For the origin 
of the PWN, it will be important to know whether there 
is any subtructure that may indicate that the nebula has 
been created as the result of a discrete number of flares, 
whether there are jets, or whether it is a more homo- 
geneous source. Deeper obscrvatons arc also important 
to obtain better spectra and images of the SNR. This 
may reveal the nature of its X-ray emission, thermal, or 
non-thermal, which will help to put better constraints 
on the age of the SNR. In our view IE 1547.0-5408 and 
the extended sources arormd it may be key objects to 
understand the origin and evolution of magnetars. 
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